home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CBAL.z / CBAL
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. ____CCCCBBBBAAAALLLL((((3333FFFF))))                                                            ____CCCCBBBBAAAALLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CBAL, SCBAL    -  EISPACK routine.  This subroutine is a complex version
  10.      of BALANCE.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ccccbbbbaaaallll((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,, aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssccccbbbbaaaallll((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh
  20.           rrrreeeeaaaallll             aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,, aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      This subroutine balances a COMPLEX matrix and isolates eigenvalues
  26.      whenever possible.
  27.  
  28.      On INPUT
  29.  
  30.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  31.      as declared in the calling program dimension statement.
  32.  
  33.      NNNN is the order of the matrix.
  34.  
  35.      AAAARRRR and AI contain the real and imaginary parts, respectively, of the
  36.      complex matrix to be balanced.  On OUTPUT
  37.  
  38.      AAAARRRR and AI contain the real and imaginary parts, respectively, of the
  39.      balanced matrix.
  40.  
  41.      LLLLOOOOWWWW and IGH are two integers such that AR(I,J) and AI(I,J) are equal to
  42.      zero if
  43.         (1) I is greater than J and
  44.         (2) J=1,...,LOW-1 or I=IGH+1,...,N.
  45.  
  46.      SSSSCCCCAAAALLLLEEEE contains information determining the
  47.         permutations and scaling factors used.  Suppose that the principal
  48.      submatrix in rows LOW through IGH has been balanced, that P(J) denotes
  49.      the index interchanged with J during the permutation step, and that the
  50.      elements of the diagonal matrix used are denoted by D(I,J).  Then
  51.  
  52.      SSSSCCCCAAAALLLLEEEE (J) = P(J),    for J = 1,...,LOW-1
  53.         = D(J,J)       J = LOW,...,IGH
  54.         = P(J)         J = IGH+1,...,N.  The order in which the interchanges
  55.      are made is N to IGH+1, then 1 to LOW-1.  Note that 1 is returned for IGH
  56.      if IGH is zero formally.  The ALGOL procedure EXC contained in CBALANCE
  57.      appears in
  58.  
  59.      CCCCBBBBAAAALLLL in line.  (Note that the ALGOL roles of identifiers K,L have been
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____CCCCBBBBAAAALLLL((((3333FFFF))))                                                            ____CCCCBBBBAAAALLLL((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      reversed.)  Questions and comments should be directed to B. S. Garbow,
  75.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.